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abstract 

We present the discovery of the optical afterglow of GRB 050505 and an optical absorption spectrum 
obtained with the Keck I 10-m telescope. The spectrum exhibits three redshifted absorption systems 
with the highest, at z = 4.2748, arising in the GRB host galaxy. The host absorption system is marked 
by a damped Lya (DLA) feature with a neutral hydrogen column density of log A^(HI) = 22.05 ± 0.10, 
higher than that of any QSO-DLA detected to date, but similar to several other recent measurements 
from GRB spectra. In addition, we detect absorption lines from both low- and high-ionization species 
from which we deduce a metallicity, Z « 0.06 Zq, with a depletion pattern that is roughly similar to 
that of the Galactic warm halo, warm disk, or disk+halo. More importantly, we detect strong absorption 
from Si II* indicating a dense environment, nn <; 10^ cm"'^, in the vicinity of the burst, with a size 
of ~ 4 pc. In addition, the C IV absorption system spans a velocity range of about lO'^ km s~^, 
which is not detected in any other absorption feature. We show that the most likely interpretation 
for this wide velocity range is absorption in the wind from the progenitor star. In this context, the 
lack of corresponding Si IV absorption indicates that the progenitor had a mass of < 25 Mq and a 
metallicity <^ 0.1 Zq, and therefore required a binary companion to eject its hydrogen envelope prior to 
the GRB explosion. Finally, by extending the GRB-DLA sample to z ~ 4.3 we show that these objects 
appear to follow a similar metallicity-redshift relation as in QSO-DLAs, but with systematically higher 
metallicities. It remains to be seen whether this trend is simply due to the higher neutral hydrogen 
columns in GRB-DLAs, or if it is a manifestation of different star formation properties in GRB-DLAs. 
Thus, GRBs hold the potential to probe all scales relevant to our understanding of the star formation 
process and its relation to metal production. 

Subject headings: gamma-rays: bursts — ISM: abundances — ISM:kinematics — stars: mass loss — 
stars: Wolf-Rayet 



1. INTRODUCTION 

The study of star formation and the associated pro- 
duction and enrichment of the interstellar medium (ISM) 
and intergalactic medium (IGM) by metals is focused on 
two observational approaches: large galaxy surveys based 
on rest-frame UV, optical, near-infrared and far-infrared 
' emission, and the use of high redshift quasars as probes of 
the ISM/IGM. The former provides a view of the star for- 
mation activity from high redshift to the present, as well as 
the possible enrichment of the IGM through galactic-scale 
winds. In the context of quasar absorption studies the 
highest column density absorbers, the so-called damped 
Lyman alpha (DLA) systems {N{m) > 2 x 10^° cm'^) 
are of particular interest. This is because of the similarity 
of their column densities to those in local luminous galax- 
ies, the fact that at these columns the hydrogen is mainly 
neutral and may form stars, and the conclusion that the 



bulk of the neutral gas in the range z — 5 is in DLAs 
(for a recent review see Wolfe et al. 2005 and references 
therein). 

The connection between luminous star-forming galaxies 
and DLAs remains an open question. This is primarily due 
to observational limitations inherent in both approaches. 
First, surveys of star- forming galaxies at z > 2 are limited 
to relatively bright object, L ^ 0.3L* (e.g., Steidel et al. 
2003). Furthermore, with the exception of the brightest 
galaxies, the continuum emission is typically too faint to 
observe absorption from interstellar gas and therefore to 
tic together the star formation and ISM properties. On the 
other hand, the use of quasar sight lines to probe the disks 
of high redshift galaxies is limited by a cross-section selec- 
tion effect to large impact parameters (^ 10 kpc). More- 
over, in the case of DLAs the small impact parameters 
severely limit the identification of the DLA counterparts 
against the bright background quasar. Thus, while some 
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DLA counterparts have been identified, primarily at z 
(e.g., Chen & Lanzetta 2003), a definitive association with 
star-forming systems or LBGs remains elusive. 

Despite the observational challenges, the various galaxy 
surveys have provided evidence for near solar mctallici- 
ties in at least some of the bright high rcdshift galaxies 
(Shapley et al. 2004; Swinbank ct al. 2004). Typical star 
formation rates in the bright galaxies are in the range of 
~ 10- 10^ Mq yr-i (Shapley et al. 2001; Forster Schreiber 
et al. 2004; Swinbank ct al. 2004). On the other hand, 
DLAs appear to be metal poor, with a typical Z ~ 0.03 
Zq (Prochaska et al. 2003), and while they exhibit evi- 
dence for star formation in a few cases, the rates typically 
appear to be lower than in the LBGs (Bunker et al. 1999). 

An alternative approach to probing intervening gas in 
galaxies and the IGM is to use the afterglows of gamma- 
ray bursts (GRBs). In the context of the relation to star 
formation and the nature of DLAs, GRBs offer several ad- 
vantages over quasar studies. First, GRBs are embedded 
in star forming galaxies with typical offsets of a few kpc 
or less (Bloom et al. 2002). They therefore not only pro- 
vide a direct link to star formation, but also probe the re- 
gions of most intense star formation, and hence production 
and dispersal of metals. Second, since the GRB afterglow 
emission fades away on a timescale of days to weeks, the 
host galaxy and any intervening DLAs can be subsequently 
studied directly (e.g., Vreeswijk et al. 2004). 

Third, and perhaps most important, since GRBs are 
likely to be located in star forming regions (e.g., molecular 
clouds) within their host galaxies, this approach provides 
the only systematic way to directly probe the small-scale 
environment and conditions of star formation at high rcd- 
shift; the probability of intersecting an individual molec- 
ular cloud in a quasar sight line is vanishingly small. In 
the same vein, GRBs can probe the circumstellar environ- 
ment of the progenitor star itself and provide a unique view 
of the mass loss history and properties of the progenitor 
(e.g., metallicity, mass, binarity). With a large statistical 
sample, this is the only way to compare the properties of 
massive stars and individual star forming regions at high 
rcdshift to those in the Milky Way and the local universe. 

Over the past several years, a few absorption spectra 
of GRB afterglows have been obtained, revealing rela- 
tively large neutral hydrogen column densities (most in 
the DLA category; e.g., Vreeswijk et al. 2004). The metal- 
licity, inferred in only a few cases, appears to be sub-solar 
{Z ~ 0.01-0.1 Zq; Vreeswijk et al. 2004; Chen et al. 2005; 
Starling et al. 2005a), but with a dust-to-gas ratio that is 
larger than that in QSO-DLAs (Savaglio ct al. 2003). In 
addition, some spectra reveal complex velocity structure, 
interpreted to arise from ordered galactic rotation (Castro 
et al. 2003), and in one case appearing to arise in the com- 
plex wind environment of the progenitor star (e.g., Moller 
et al. 2002). These results already suggest that GRBs 
probe environments that are likely missed in the quasar 
surveys. 

With the advent of the Swift satellite, we can now start 
to utilize GRBs as probes of the high redshift universe 

in a systematic manner. Swiff s sensitivity and ability to 
rapidly and accurately localize a large number of GRBs 



has resulted in a redshift distribution spanning nearly uni- 
formly from z ~ 0.5 to ^ 6 (Bergcr ct al. 2005; Jakobsson 
et al. 2005b; Kawai et al. 2005; Haislip et al. 2005). This 
sample is therefore well-matched to the star formation his- 
tory of the universe, and over time will allow us to address 
the redshift evolution of star forming environments and 
perhaps individual massive stars. 

Along these lines, we present here an absorption spec- 
trum of GRB 050505, which reveals a DLA with a col- 
umn density logiV(HI) = 22.05 ± 0.10 at a redshift of 
z = 4.2748. This system is currently the highest redshift 
GRB host for which detailed information is available. The 
spectrum probes not only the interstellar medium of the 
host galaxy, but also provides information on the local 
environment of the burst, likely including the wind of the 
progenitor star from which we are able to draw conclusions 
about the nature of the star that exploded. 

2. OBSERVATIONS 

GRB 050505 was detected by Swift on 2005 May 5.974 
UT. The duration and fluence of the burst are 60 s and 
(4.1 ± 0.4) X 10~^ erg cm~^, respectively (HuUinger et al. 
2005). Observations with the Swift X-ray telescope (XRT) 
started 47 min after the burst, and revealed an uncata- 
logued source at a=09'^27'"03.2% 5=-F30°16'21.5" (J2000) 
with an uncertainty of about 6" and a flux of~2xl0~^^ 
erg cm~^ s""'^ (Kennea et al. 2005) . No object was de- 
tected by the Swift UV/optical telescope (UVOT) in the 
first 8 hours to a limit of F > 20.35 mag and B > 21.04 
mag, at a mean time of 2.49 and 2.59 hr, respectively 
(Rosen et al. 2005). 

We initiated observations of GRB 050505 with the Low 
Resolution Imaging Spectrometer (LRIS) mounted on the 
Keck I 10-m telescope about 6.4 hr after the burst. We 
obtained simultaneous g- and /-band observations and de- 
tected an object close to the center of the XRT error cir- 
cle at a=09'^27'"03.3^ (5=+30°16'23.7" (J2000), with a 
brightness / = 20.51 ± 0.05 mag and g ~ 23.67 ± 0.12 mag 
(Figure 1). The spectral slope between the two bands, 
oc jy-4-9±o.3^ jg sharp for host galaxy extinction, 
and instead suggests a redshift 2; ~ 4 — 5.5. Contempora- 
neous observations with UKIRT revealed a near-IR coun- 
terpart with a brightness of if = 18.1 ± 0.2 mag (Rol et al. 
2005). 

Following the identification of the afterglow we used 
LRIS to obtain two 900-s spectra with a 1" wide slit. 
The spectra were reduced using standard IRAF routines, 
while rectification and sky subtraction were performed us- 
ing the method of Kelson (2003). Wavelength calibra- 
tion was performed using HgArNcZnCd arc lamps and 
air-to-vacuum and heliocentric corrections were applied. 
The resulting dispersion scales are 1.86 A/pix for the red 
side and 0.61 A/pix for the blue side, with an rms wave- 
length uncertainty of about O.2A. Finally, flux calibra- 
tion was performed using the spcctrophotometric standard 
BD+284211, with a correction for the small amount of 
Galactic extinction, E{B — V) = 0.021 mag (Schlegel et al. 
1998). Convolution of the spectrum with the / filter band- 
pass results in a flux density of 15.8 ± 1.5 fiJy (20.5 ± 0.1 
mag), in very good agreement with the measured 7-band 
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magnitude. This indicates that slit losses were minimal. 

We identify three redshift systems in the extracted spec- 
trum, zi = 4.2748 ± 0.0008,' Z2 = 2.2650 ± 0.0008, and 
zz = 1.6948 ± 0.0003. We associate the redshift system z-^ 
with the host galaxy of GRB 050505 given the nature of the 
absorption lines (see §3) and the lack of Lya forest absorp- 
tion redward of the damped Lya feature at the redshift zi . 
At this redshift, using the standard cosmology {Hq = 71 
km s~^ Mpc~^, flrn = 0.27 and Oa = 0.73), the isotropic- 
equivalcnt 7-ray energy release is ii^-,,iso ~ 8.9 x 10"'''^ erg, 
while the X-ray luminosity extrapolated to t = 10 hr (as- 
suming the typical Fi, oc t~^-^) is Lx,iso ~ 2.3 x 10"*^ erg 
s~^ (Nousek et al. 2005). Both values are at the high end 
of the distribution for previous GRBs (Berger et al. 2003; 
Bloom et al. 2003). 

3. THE HOST GALAXY ABSORPTION SYSTEM 

The spectrum of GRB 050505 is dominated by a broad 
absorption feature centered on an observed wavelength of 
about 6400A, which we identify as Lya. A Voigt profile fit 
to the Lya absorption feature results in a column density 
of neutral hydrogen, log jV(HI) = 22.05 ± 0.10, indicating 
that the absorber is a DLA, with one of the highest column 
densities measured to date in either QSO or GRB sight 
lines (Curran et al. 2002; Vreeswijk et al. 2004). We note 
that the actual column density may be somewhat higher 
since, unlike in the case of quasars, the GRB is embedded 
in the host galaxy (Vreeswijk et al. 2004). The average 
flux decrement between Lya and Ly/3 is about 0.55, which 
is in good agreement with the values measured in quasar 
sight lines at the same redshift (Cristiani et al. 1993). 

Blueward of the damped Lya feature we detect absorp- 
tion from the Lya forest, likely mixed with metal absorp- 
tion features of oxygen, silicon, nitrogen, sulfur, carbon, 
and iron. Due to the low resolution of our spectrum we 
cannot disentangle the contribution of the metals from the 
forest. We also detect Ly/3 absorption, and at a wavelength 
of about 4880A we detect the Lyman limit. 

Redward of the Lya absorption feature we detect a wide 
range of metal lines. The identifications, observed wave- 
lengths, and equivalent widths (EWs) of these lines are 
listed in Table 1. We note that in some cases these lines 
are blended with absorption features from the lower red- 
shift system at Z2 = 2.265. Most of the lines detected in 
the spectrum appear to be saturated. We can obtain a 
lower limit on the column density assuming the optically 
thin case 
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where / is the oscillator strength, W\ is the equivalent 
width, and A is the rest wavelength. The resulting lower 
limits for each transition are listed in Table 1. 

Since most of the lines are saturated and therefore do 
not lie on the linear part of the curve of growth (COG) 
we construct a joint COG for all the detected transitions 
which do not suffer from strong blending. In this case we 
use the standard formulation (Spitzer 1978) 
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We assume that all species share the same value of the 
Doppler parameter, b, and we fit iteratively for b and the 
column density of each species (e.g., Savaglio et al. 2003). 
The resulting best-fit COG is shown in Figure 3, and the 
column densities are listed in Table 2. We find that lines 
of C I and Ni II are well-described by the linear part of the 
COG, lines of S II are mildly saturated, and most other 
species lie on the flat portion of the COG. We stress that 
given the assumption of a single b value, as well as the low 
resolution of our spectrum, the derived column densities 
of the strongly saturated lines should still be considered 
as rough lower limits. 

With this caveat in mind we derive the following abun- 
dances from the COG analysis. The column density of 
S II is logiV « 16.1, leading to an abundance relative to 
the solar value of [S/H] » —1.2; here we use the compila- 
tion of solar abundances of Asplund et al. (2005). Sulfur 
is a non-refractory element and its gas-phase abundance 
therefore closely matches the gas metallicity. The derived 
value is roughly similar to those of both QSO- and GRB- 
DLAs. The Fe II column density of log A'' w 15.5 indi- 
cates an abundance of [Fe/H] « —2.0, which is similar to 
the mean value for QSO-DLAs at a similar redshift range 
(e.g., Prochaska et al. 2003). The ratio of sulfur to iron, 
[S/Fe] w 0.8, however, is at the high end of values mea- 
sured for QSO-DLAs, which at [Fe/H] w —2.0 range from 
about 0.1 to 1 (Lopez & Ellison 2003). Since iron is a 
strongly depleted element, the large value of [S/Fe] sug- 
gests a high dust content, similar to the inference made in 
several other GRB-DLAs (Savaglio et al. 2003; Savaglio & 
Fall 2004; Watson et al. 2005). 

Inspecting the abundances of other ions we find that 
[Si/H] « —1.6, when we sum the contributions of Si II and 
Si IV. The ratio compared to iron is [Si/Fe] « 0.4, which 
is in excellent agreement with the median value of about 
0.4 for QSO-DLAs (Wolfe et al. 2005). 

Using the abundances of sulfur, silicon, and iron we pro- 
vide a comparison to the Milky Way depletion patterns in 
Figure 4. We consider the four typical patterns of warm 
halo (WH), warm disk-hhalo (WDH), warm disk (WD), 
and cool disk (CD) clouds (Savage & Sembach 1996), and 
follow the method of Savaglio et al. (2003) to determine 
the metallicity and dust-to-gas ratio in each case. We find 
that the first three depletion patterns provide an adequate 
representation of the data with a metallicity Z « 0.06 Zq 
and dust-to-gas ratios of 1.1 (WH), 0.95 (WDH), and 0.85 
(WD). The cool disk depletion pattern does not provide 
an adequate fit since it underestimates the silicon abun- 
dance and over-estimates the iron abundance. Given the 
inferred metallicities and dust-to-gas ratios we calculate a 
rest-frame F-band extinction of about 0.3 ± 0.1 mag. 
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3.1. A CW Outflow 

The most prominent metal absorption feature in the 

spectrum of GRB 050505, with a total rest- frame equiva- 
lent width of about 11.6A, is a blend of C IVAA 1548, 1550. 
The overall velocity spread of this feature is about 950 km 
s~^, extending blueward of the systemic redshift of the 
host galaxy, as defined by other metal absorption features 
(Figure 5). The wide velocity spread is not observed in 
any of the other low- or high-ionization features; for ex- 
ample the limit on blue-shifted Si IV is log A'' < 13.5, at 
least a factor of thirty below the C IV column density. It 
is not clear from the low resolution spectrum if the veloc- 
ity structiu'c is due to a set of discrete absorbers or to a 
relatively uniform distribution of C IV. 

Regardless of this distinction, the observed velocity 
range is much larger than the typical values associated 
with the rotation or velocity dispersion of individual galax- 
ies, including those observed in past GRB absorption spec- 
tra (e.g., GRB 000926 with v « 170 km s"^; Castro 
et al. 2003). Similarly, blueshifted absorption by inter- 
stellar low- and high-ionization gas has been detected in 
the stacked spectra of Lyman break galaxies (LBGs) , but 
with a typical velocity shift relative to the nebular emis- 
sion lines of only ~ —150 km s~^ (Adelberger et al. 2003). 
On the other hand, Frye et al. (2002) find evidence for 
outflows in z ^ 4 lensed galaxies which approach 800 km 
s~^, but these outflows are evident in low- ionization lines 
of O I, Si II, and C II. In local bright starburst galax- 
ies outflow velocities are typically < 300 km s""'^ (Heck- 
man et al. 2000). Therefore, unless the wind from the 
host galaxy of GRB 050505 is unusually fast, we consider 
this interpretation unlikely. This would particularly be the 
case if the host of GRB 050505 is similar to those of other 
GRBs, which tend to have relatively low masses (Chris- 
tensen et al. 2004) and would therefore have escape veloc- 
ities much smaller than 10^ km s~^. 

The C IV absorption may be alternatively related to 
an extended large-scale structure along the line of sight. 
Adelberger et al. (2003) find a strong correlation between 
the location of C IV absorbers in quasar sight lines and 
the location of LBGs. In particular, for a C IV column 
density of > 10^** cm^^ (as observed in the spectrum of 
GRB 050505) about 80% of the absorbers lie within 600 km 
s~^ and A9 = 35" of an LBG. While this velocity range is 
similar to the one observed here, it has been argued that 
the actual outflow velocities are ^ 400 km s~^ (Songaila 
2005), and that the extent of the C IV bubbles may be a 
reflection of pre-galactic enrichment at z ~ 10 by slower 
winds from dwarf galaxies (Porciani & Madau 2005). Since 
here we directly detect an outflow with v ^ 10^ km s^^ it 
is unlikely to be related to the C IV bubbles discussed in 
Adelberger et al. (2003) and Songafla (2005). 

Another possibility in the context of an extended struc- 
ture is that the C IV absorption arises from an overlapping 
structure of galactic lialos, perhaps stretched out along a 
dark matter filament. This hypothesis is difficult to assess 
with a single absorption spectrum. However, if such struc- 
tures are ubiquitous they should be observed in a large 
fraction of GRB absorption spectra, as well as in quasar 
absorption spectra. Studies of C IV systems in quasar 



spectra suggest that these systems are strongly clustered 
for velocity shifts of ^ 200 km s^^ and are essentially un- 
correlated for v ^ 500 km s^^ (Ranch et al. 1996; Pichon 
et al. 2003). Thus, we conclude that the observed C IV ve- 
locity structure is not due to a galactic-scale phenomenon. 

In the context of a massive star progenitor, an attrac- 
tive alternative is that the C IV velocity structure is the 
signature of a fast wind. A similar inference was made in 
the case of GRB 021004, which exhibited a complex ve- 
locity structure extending from about 140 to 3000 km s^^ 
(Moller et al. 2002; Mirabal et al. 2003; Schacfer et al. 
2003; Fiore et al. 2005; Starhng et al. 2005b). The pres- 
ence of such a fast outflow led to the general conclusion 
that the progenitor was a Wolf-Rayet star. However, un- 
like in the case presented here, the velocity structure in 
the spectrum of GRB 021004 was evident in both low- and 
high-ionization lines, including H I. This raised the diffi- 
culty of explaining the presence of both types of species in 
the highly- ionized burst environment, as well as the pres- 
ence of a large column of hydrogen in a Wolf-Rayet wind. 

While these difficulties prevented a definitive explana- 
tion of the high-velocity absorbers, two main possibilities 
have been discussed. Mirabal et al. (2003) argue that to 
reach the observed expansion velocities, radiative accelera- 
tion of the pre-existing fragmented shell nebula is required. 
In this model the interaction of the fast Wolf-Rayet wind 
with the slower winds of earlier mass loss episodes results 
in a fragmented structure due to Rayleigh- Taylor instabil- 
ities. The radiation pressure produced by the GRB and 
afterglow emission then accelerates the fragmented shell 
nebula to velocities of hundreds to thousands of km s^^. 

van Marie et al. (2005), on the other hand, produce hy- 
drodynamic numerical simulations of the interactions of 
various mass-loss phases (fast main sequence wind, slow 
red supergiant wind, and a fast wolf-Rayet wind), and 
conclude that these interactions lead to a clumpy struc- 
ture with a velocity range in excess of 2000 km s~^. The 
fragmented and clumpy structure of the wind is a common 
theme in both scenarios, and is used to explain the similar 
kinematic structure in both the low- and high-ionization 
lines. 

In the case of GRB 050505 the velocity structure is 
only apparent in C IV, suggesting that a more uniform 
wind structure is allowed. Still, the effect of the ioniz- 
ing radiation of the GRB has to be taken into account. 
Since the recombination tinicscalc is expected to be sig- 
nificantly longer than the lifetime of the GRB (unless 
the density is ^ 10^ cm~^; Perna & Loeb 1998), we 
can estimate the size of the ionized region roughly as 

r ~ ^/ -Bioncrciv / 47reion,civ ~ few x 10^^ cm; here E'ion 
is the total energy in photons that can ionize C IV, 
cciv = 1-068 X 10^^® cm^ is the photoionization cross- 
section, and Cion.civ = 64.5 eV is the ionization threshold 
energy. Clearly, a detailed photoionization calculation is 
required, but we conclude that outside of a few parsec C IV 
ions should remain intact. We stress that the minimum ra- 
dius for survival of Si IV ions is a factor of about 3.5 times 
smaller than that of C IV, so that the non-detection of 
outfiowing Si IV cannot be simply attributed to photoion- 
ization. 
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Thus, the observations require a fast-moving 10^ km 
s~^) wind, enriched in carbon and deficient in silicon, and 
with a radius of at least a few parsec. These requirements 
can be naturally satisfied in the case of a Wolf-Rayet wind 
for which typical velocities can easily exceed 10^ km 
and extend out to a distance of ^ 20 pc (Abbott 1978; van 
Marie et al. 2005). 

The lack of silicon absorption, with a column density 
that is at least a factor of thirty lower than that of C IV, 
provides interesting constraints on the nature of the pro- 
genitor. Calculations of theoretical profiles of C IV and 
Si IV in the winds from a range of massive stars suggest 
that while strong C IV absorption is nearly always present, 
the presence of Si IV lines depends sensitively on a star's 
mass and metallicity (Leitherer & Lamcrs 1991). Thus, 
the wind from progenitors with a low mass and metallicity 
will have negligible Si IV absorption. For example, in the 
case of a wind from a 25 Mq star at the end of core hydro- 
gen burning the ratio of Si IV to C IV matches the observed 
limit oiEWiSilY)/EW{ClY) < 0.1 for Z < 0.1 Zq, while 
for a 60 Mq star the theoretical ratio is ~ 0.5 even for 
Z = 0.03 Zq. These considerations indicate that the pro- 
genitor of GRB 050505 was most likely a carbon-rich (WC) 
Wolf-Rayet star with a relatively low mass 25 Mq) and 
metallicity ^ 0.1 Zq). 

However, radiatively-driven winds follow a mass loss 
rate metallicity relation, M oc ^0-5~o.7 (Maeder & Meynet 
1988; Leitherer & Lamers 1991) and therefore single stars 
reach the Wolf-Rayet phase at a higher mass when their 
metallicity is low; in the Galaxy MwR.min ~ 25 Mq , while 
in the SMC AfwR.min « 45 M© (Maeder 1998; Maeder 
& Meynet 2000). With a metaUicity of ^ 0.1 Zq and a 
mass of only ^ 25 Mq, it appears that the progenitor of 
GRB 050505 was not sufficiently massive to have become a 
Wolf-Rayet star. It is therefore likely that the progenitor 
required a companion star in order to eject the hydrogen 
envelope. 

3.2. Sill* 

Another interesting feature of the spectrum of 
GRB 050505 is the detection of strong fine-structure ab- 
sorption from Si II*. Since these transitions have not been 
convincingly detected in QSO-DLAs (e.g., Howk et al. 
2005) and since they are most likely excited by collisions 
with electrons (for T < 10^ K; Silva & Viegas 2002), we 
conclude that the Si II* absorber is coincident with the 
local, high density environment of the GRB. In Figure 6 
we show the three detected Si II and Si II* transitions. 
We note that the Si IIA1304 and Si II*A1309 pair is lo- 
cated in the atmospheric B-band, therefore precluding a 
clear measurement of the equivalent widths. In addition, 
due to the low spectral resolution, the Si IIA1260 transi- 
tion is blended with S IIA1259 and Fe IIA1260, although 
given the large oscillator strength of the Si II transition, 
it likely dominates the measured equivalent width. Our 
estimate of the column density is thus based primarily on 
Si IIA1527 and on Si II*AA1264, 1533 (Figure 6). 

The absorption lines of both species are saturated and 
the ratio of column densities is therefore somewhat un- 
certain, although it is clear that the column of Si II* is 



fairly large; using the weakest, unblended line of Si II* 
under the assumption of the optically thin case we find 
log A^ (Sill*) 14.7. This value is somewhat higher than 
those found in two previous GRB spectra: log A/' « 14.2 
(GRB 030323; Vreeswijk et al. 2004) and logA^ « 14.3 
(GRB 020813; Savagho & Fall 2004), but the rough sim- 
ilarity may reflect comparable physical conditions in the 
local environments of GRBs. 

Using the COG analysis we find that the ratio of column 
densities is 

^«1^«02 (5) 

With this value, and using the calculations of Silva & Vie- 
gas (2002), we estimate that for a temperature of 10^ K 
(i.e., a neutral medium), the implied volume density of H I 
is nui ^ 10^ cm~^ if the electron fraction is = O.lnni- 
If the electron fraction is rie < lO^^rifji, as may be ex- 
pected for a neutral medium in which the electrons pri- 
marily come from low-ionization species with an abun- 
dance relative to hydrogen of ~ 10~^ (Silva & Viegas 
2002; Vreeswijk et al. 2004), the density is jihi 
cm~^. In this framework, the size of the absorbing region 
is £ ~ 10^^ °^/10^ ~ 3.5 pc assuming that the hydrogen 
is co-located with the Si II* absorber. The inferred mass 
of the absorbing region is M = mpN(Hl)£'^ ~ 10'^ Mq. 
Naturally, if only a fraction of the H I column is associ- 
ated with the Si II* absorber, then the derived mass and 
size are in fact upper limits. The inferred compact size 
of the absorber explains the non-detection of Si II* in the 
QSO-DLAs, since the probability of intersecting a region 
of only a few parsec in size is vanishingly small. 

We note that there are two difficulties with this interpre- 
tation if we insist that the compact Si II* absorber has to 
be placed around the GRB. First, in order to survive the 
burst of ionizing radiation from the GRB (which affects 
the environment on a scale of ~ 10 pc) the Si II absorber 
has to be shielded. Second, as discussed in the previous 
section, the wind from the progenitor star appears to be 
deficient in silicon, and another source for the Si II may 
be required. These arguments suggest that the Si II is 
located in the vicinity of the GRB, but is probably not 
directly associated with the progenitor star itself, and is 
instead a signature of the dense environment in the star 
forming region. 

3.3. Limits on the H2 Column Density 

Absorption lines of molecular hydrogen may be detected 
in the spectrum of GRB 050505 if the burst occurred in a 
dense molecular cloud. Unfortunately, at the low resolu- 
tion of our spectrum it is essentially impossible to discern 
H2 absorption lines from the Lya forest and transitions 
due to metals. We note that in general the molecular 
fraction in DLAs is <J 10~^, and only a small number 
of DLAs exhibit detectable fractions of H2 (Ledoux et al. 
2003). The explanation for this trend is a combination of a 
low dust content and fairly intense far-UV radiation field, 
which combine to reduce the formation rate of H2 and 
increase the photodissociation rate, respectively (Wolfe 
et al. 2005). However, if GRBs tend to occur in molec- 
ular clouds, GRB-DLAs may exhibit preferentially higher 
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fractions of molecular hydrogen, particularly in light of the 
fairly largo dust-to-gas ratios in comparison to QSO-DLAs 
(Savaglio et al. 2003). 

On the other hand, as discussed in detail by Draine & 
Hao (2002), the optical/UV radiation from the GRB itself 
is expected to photoionizc and dissociate H2 molecules to 
a distance of a few parscc. If the cloud hosting the GRB 
progenitor is smaller than this size, or if the burst is lo- 
cated near the edge of the cloud we expect that no H2 
absorption features will be detected even if the molecular 
fraction was initially high. 

If the cloud is sufficiently large, the UV radiation 
field is expected to also excite H2 into vibrational levels 
which will produce observable absorption features at rest- 
wavelengths of 1110 - 1705A (Draine & Hao 2002). The 
strongest features are at A = 1254. SA and A = 1277. SA. 
While we detect an absorption feature coincident with the 
I277.3A line, we interpret this feature as C IA1277.2 whose 
equivalent width is in good agreement with other C I lines. 
The non-detection of vibrational H2 absorption, under the 
assumption that the cloud size is sufiiciently large, places 
a limit of log A''(H|) ^ 18 cm'^ (Draine & Hao 2002), and 
hence /(H2) ^ 10"*. 

3.4. Limits on Lya Emission 

We do not detect Lya in emission with a 3a limit of 
F < 2.1 X 10"^'^ erg cm-^ s'^. At the redshift of the 
burst this limit translates to an upper limit on the line 
luminosity of _L < 3.9 x 10^^ erg s^^. Using the conversion 
of Kennicutt (1998) and assuming case B recombination, 
the limit on the star formation rate implied by the non- 
detection of Lya emission is SFR = 9.1 x 10~^^IyLya < 3.5 
Mq yr-i. 

The limit on the inferred star formation rate is about 
a factor of three higher than the value inferred for 
GRB 030323 at ^ = 3.372 (Vreeswijk et al. 2004), but 
it is significantly lower than the value measured for 
GRB 021004, sa 11 Mq yr-\ based on Lya omission (e.g., 
Fynbo et al. 2005). In fact, nearly every GRB host in 
which Lya emission can be detected has been detected 
(Jakobsson et al. 2005a), while the fraction of LBGs with 
Lya equivalent widths similar to those in GRB hosts is 
only ~ 1/3 (Shapley et al. 2003). 

We note that the lack of Lya emission may still accom- 
modate a star formation rate higher than our formal limit 
since the Lya photons are easily destroyed by resonant 
scattering in the presence of dust. Thus, even for a fairly 
low dust content it is possible that the star formation rate 
is in fact higher than the limit provided above. A bet- 
ter assessment of the star formation rate requires near-IR 
spectroscopy of the host galaxy for detection of Ha at 3.46 
fim and/or [O II]A3727 at 1.97 /xm. 

4. INTERVENING SYSTEMS 

We detect two intervening systems at Z2 = 2.265 and 
Z3 = 1.695. The latter is detected only in Mg II with a 
rest-frame equivalent width of 1.98A (A2796.35) and 0.94A 

(A2803.53). These values match the expected 2 : 1 ratio 

®An absorption spectrum of GRB 050904 at z = 6.29 has been ob1 
not been published so far. 



in the optically thin case, suggesting that the lines are not 
saturated. The derived column density is about 5 x lO^'' 
cm~^. A comparison to the sample of Steidel & Sargent 
(1992) reveals that 10% of quasar sight lines reveal sys- 
tems with similarly large equivalent widths. However, un- 
like the observed doublet ratio of about two measured here, 
the range of doublet ratios for the quasar Mg II systems 
with EW ^ 2 A is typically 1 — 1.3, indicative of strong sat- 
uration. We finally note that for a magnesium abundance 
of about 10% the solar value, the inferred column density 
in the Z3 system suggests an associated neutral hydrogen 
column density of logiV(HI) ~ 18.5 cm~^. 

The Z2 redshift system exhibits a wider range of ab- 
sorption features, but unfortunately, some of these lines 
are blended with features from the host galaxy. Still, sev- 
eral Fe II lines, as well as the Mg II doublet, and a sin- 
gle Mn II line are not blended, allowing us to estimate 
the column density of some metals in the Z2 system. A 
COG fit to the three Fe II lines indicates 6 « 60 km s~^ 
and log AT (Fell) w 15.0. The Mg II ratio of about 0.85 
indicates strong saturation, and a lower limit on the col- 
umn density of logA''(MgII) « 14.6, assuming the same 
Doppler parameter as for the Fe II lines. Finally, assum- 
ing that the single Mn II line is optically thin we esti- 
mate logiV(Mnll) w 13.7. Thus, the relative abimdance 
is [Fe/Mn] ^ -0.7, while in DLAs it is typically ~ 0. 
This suggests that we may be underestimating the Fe II 
by about 0.7 dex. 

We also note that both iron and manganese are strongly 
depleted in the Galaxy, while elements such as zinc are rel- 
atively undepleted (Savage & Sembach 1996). In our case, 
the Zn IIA2062.66 line is blended with C I from the host 
galaxy system. However, using the measured equivalent 
width as an upper limit and assuming the optically thin 
case, we find logAr(Znll) < 13.4. Thus, [Mn/Zn] > -0.5 
and [Fe/Zn] > —1.2 are in agreement with the typical 
values measured for DLAs of about —0.5 and —0.4, re- 
spectivcily, particularly if the abundance of Fe II is in fact 
underestimated as suggested above. 

5. DISCUSSION AND CONCLUSIONS 

The absorption spectrum of GRB 050505 exhibits a 
DLA with a column density of log A^ (HI) « 22.05 associ- 
ated with the host galaxy. The measured column is higher 
than that of any QSO-DLA studied to date, but is similar 
to those measured for several GRB-DLAs, albeit at the 
highest redshift to date^. In Figure 7 we provide a sum- 
mary of the column densities measured from both quasars 
and GRBs to date. As already noted in Vreeswijk et al. 
(2004), these distributions differ significantly, presumably 
as a result of the smaller impact parameters in the case of 
GRB-DLAs and the possibility that they probe individual 
dense star forming regions. The metallicity we infer for 
GRB 050505, Z w 0.06 Zq, is within the range of values 
measured for QSO-DLAs at a similar redshift (Prochaska 
et al. 2003), but is a factor of several higher than the mean 
metallicity of about 0.01 Zq. 

This, along with the metallicities inferred for five addi- 

by Kawai et al. (2005), but the column density and metallicity have 



7 



tional GRB-DLAs, is shown in Figure 8. In the top panel 
we plot the metallicity as determined from various indica- 
tors as a function of redshift in comparison to QSO-DLAs. 
For the QSO sample, Prochaska et al. (2003) noted an ap- 
parent trend of increased metallicity with lower redshift. 
A similar trend appears in the admittedly small GRB- 
DLA sample which spans from z w 2 to 4.3. However, 
it is also clear that the GRB-DLAs have systematically 
higher metallicities compared to QSO-DLAs. In the bot- 
tom panel of Figure 8 we plot the metallicity as a function 
of hydrogen column density. No clear trend is evident, per- 
haps suggesting that the higher metallicities in GRB-DLAs 
may not be simply a function of the column density, but 
may reflect a difference in the star formation properties, 
rates, or history of GRB-selected host galaxies compared 
to QSO-DLAs. 

We highlight two aspects of the spectrum, which provide 
insight into scales and conditions that are not typically 
probed in quasar sight lines: A C IV outflow that is most 
likely associated with the wind of the progenitor star, and 
strong absorption from fine-structure Si II* that is associ- 
ated with a dense region in the vicinity of the burst. The 
latter suggests that the GRB progenitor lived in an envi- 
ronment typical of a molecular cloud. The C IV outflow. 



taken in association with the lack of corresponding Si IV 
absorption suggests that the progenitor star had a rela- 
tively low mass and metallicity, M ~ 25 and Z <^^.\ 
Zq. The combination of these properties suggests that the 
progenitor was part of a binary system, with the compan- 
ion responsible for stripping the hydrogen envelope. 

The continued detection of high-velocity features in 
GRB absorption spectra may thus hold the only key to 
a detailed understanding of the distribution of progenitor 
properties (mass, metallicity, binarity), as well as the abil- 
ity to probe the properties of individual molecular clouds 
at high redshift. Coupled with information on the galac- 
tic scale properties of the ISM, GRBs can probe all the 
scales relevant to our understanding of the star formation 
process and its relation to metal production, as well as the 
relation between DLAs and star forming galaxies. 
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Line Identification 
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Note. — Absorption features identified in the spectrum of GRB 050505. We do not include metal lines 
blueward of the Lya absorption since at the low resolution of our spectrum these lines are blended with 
Lya forest features. The columns are (left to right): (i) Observed wavelength, (ii) line identification, 
(iii) oscillator strength (Prochaska et al. 2003), (iv) redshift of the line, (v) rest- frame equivalent width, 
and (vi) logarithm of the column density assuming the optically-thin case (Equation 1); in most cases 
this is a lower limit since the lines are genearlly saturated. In the case of blended features we derive 
the column density by assuming that the total equivalent width is due to the highest redshift feature 
(typically the host galaxy). " Absorption features are located in the atmospheric B band; Absorption 
features are located in the atmospheric A band. 



Table 2 

Column Densities of Ions in zi 



Ion 


log AT 
(cm-3) 


[X/H] 


HI 


22.05 




CI 


14.6 


-3.8 


CII 


17.0 


-1.4 


C IV 


17.1 


-1.3 


O I 


16.5 


-2.2 


Al II 


16.3 


-0.3 


Si II 


15.7 


-1.8 


Si II* 


15.1 


-2.4 


Si IV 


15.6 


-2.0 


SII 


16.1 


-1.2 


Fell 


15.5 


-2.0 


Ni II 


14.6 


-1.7 



Note. — Ionic column densities and abun- 
dances relative to the Solar values as de- 
rived from the curve-of-growth analysis 
(Figure 3). These values arc lower limits 
in most cases given the low spectral resolu- 
tion and the location of most lines on the 
flat portion of the COG. Thus, the metallic- 
ity as inferred from the sulfur abundance is 
> 0.06 Zq. 




Fig. 1. — Discovery image of the optical afterglow of GRB 050505 obtained with LRIS on the Keck I 10-m telescope. The 
inset is a combined g + I flux-calibrated color image clearly showing the red color of the afterglow which is due to the 
damped Lya absorption, the Lya forest, and the Lyman hmit absorption at A ^ 4900A (see Figure 2). 
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Fig. 2. — Absorption spectrum of GRB 050505 obtained with LRIS on the Keck I lO-m telescope. Observational details 
are given in §2. Metal absorption features from all three systems are shown, including lines that are blends from both zi 
and Z2- We note that absorption features blueward of the damped Lyman alpha feature are strongly blended with the 
Lyman alpha forest at the resolution of our spectrum and their position is shown only for completeness. The inset shows 
a zoom-in of the Lyman alpha absorption. The solid line is the best fit with log N (HI) = 22.05, and the dashed lines 
designate the la uncertainty of 0.1 dex. 
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Fig. 3. Curve of growth (COG) for the host galaxy system of GRB 050505. We constructed the COG by itcratively 
fitting for the column densities of individual ions and the Doppler parameter, b, which we assumed to have a single value. 
Transitions in the linear part of the COG lead to well-determined columns (modulo the low resolution of our spectrum). 
Transitions on the flat portion of the COG are sensitive to the value of b and should be considered as lower limits. 




B Si Mg Mn Cr Fe Ni 



Fig. 4. — Depletion pattern for the host galaxy of GRB 050505 using the inferred columns of sulfur, silicon, and iron 
relative to hydrogen. The lines represent the depletion patterns expected for various phases of the interstellar medium 
based on observations in the Milky Way (Savage & Sembach 1996), with the dust-to-gas (k) and metallicity left as free 
parameters. We find that a cool disk interpretation is unlikely. 
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Fig. 5. Absorption features of C IV and Si IV. The C IV doublet exhibits a velocity spread of about 1000 km s~^, 
which is not detected in Si IV. The gray line shows the expected absorption profile from two systems at zia = 4.2741 and 
ziB = 4.2572 with C IV column densities of 1 x 10^^ and 2 x 10^^ cm~^, respectively, and 6 = 90 km s~^ (see Figure 3). 
This scenario provides an adequate fit to the data, and indicates that the velocity structure is most likely due to a fast 
wind from the progenitor star. 
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Fig. 6.— Absorption lines of Si II and the fine structure Si II*. The hues of both species appear to be saturated so the 
ratio of column densities is only roughly determined to be ~ 0.2. We note that Si IIA1304 and Si II*A1309 are located in 
the atmospheric B band and the derived equivalent widths are highly uncertain. 
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Fig. 7. — Fractional distribution of neutral hydrogen column densities from quasar (Curran et al. 2002) and GRB 
observations (Jakobsson et al. 2004; Vreeswijk et al. 2004; Chen et al. 2005; Prochaska et al. 2005b, a; Starling et al. 2005a; 
Vreeswijk et al. 2005; Watson et al. 2005). The GRB-DLAs exhibit higher column densities with about a third exceeding 
the values measured in any QSO-DLA. 
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Fig. 8. Top: Mctallicity as a function of redshift for QSO-DLAs (gray; Prochaska et al. 2003) and GRB-DLAs (black; 
Savaglio et al. 2003; Vreeswijk et al. 2004; Chen et al. 2005; Starling et al. 2005a; Vreeswijk ct al. 2005; Watson et al. 
2005); GRB 050505 is marked with a cross. Bottom: Metallicity as a function of neutral hydrogen column density. The 
trend of increased metallicity with lower redshift, noted for QSO-DLAs by Prochaska et al. (2003), is also apparent in the 
GRB-DLA sample, albeit at an overall higher metallicity. This suggests that GRBs select more metal-enriched systems, 
but it is not clear at the present if this is related to the higher column densities. 



